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Predicted by hierarchical structure formation, Milky Way-type galaxies are anticipated to host numerous dark 
matter subhalos with masses between 10^" and a cut-ofT of 1O~^M0, or even less. In self-annihilating dark 
matter scenarios, these objects could be visible in the 7-ray band as faint and non-variable sources without 
astrophysical counterpart. In accordance with realistic subhalo models and current observational constraints 
on self-annihilating dark matter, we predict that about one massive Galactic subhalo may have already been 
detected in the 11-months catalog of Fermi-hAT data (IFGL). Selection cuts applied to the IFGL reveal twelve 
possible candidates, and in-depth studies of the most promising object, IFGL J0030. 7-1-0724, are presented. In a 
dedicated X-ray follow-up observation with the Swift XRT, seven point-like X-ray sources have been discovered. 
Within the positional uncertainty derived from the 24-months Fermi-hA'T data, we consider the unidentified 
radio source NVSS J003119-(-072456, coinciding with one of the discovered Swift sources, as the most promising 
counterpart candidate for IFGL J0030. 7-1-0724. The broad-band spectral energy distribution is consistent with 
a high-energy-peaked blazar. However, flux and extent of the Fermi source may also be compatible with a dark 
matter subhalo. A discrimination between the two scenarios requires further multi-wavelength observations. 
Strategies for identifying 7-ray sources associated with self-annihilating DM subhalos are discussed. 



1. Introduction 

Unraveling the nature of dark matter (DM) is part 
of the major tasks in modern astro- and particle 
physics. Various independent astrophysical observa- 
tions indicate a non-baryonic form of cold dark matter 
(CDM) to prevail over the baryonic content of the Uni- 
verse, e.g., [HQ]- Provided by theories extending the 
standard model of particle physics (SM), such as su- 
persymmetry and universal extradimensions, a class of 
promising CDM candidates are stable weakly interact- 
ing massive particles (WIMPs). By self-annihilation 
in SM final states (heavy quarks, gauge bosons, or 
leptons), WIMPs can produce detectable signatures 
such as photons, antimatter, and leptons, arising from 
hadronisation and subsequent decay of the final an- 
nihilation states. Apart from (loop-suppressed) line 
contributions, these processes result in continuous 
photon spectra, following a hard power law (index 
r < 1.5) cutting off exponentially to the WIMP mass. 

Within the framework of hierarchical structure for- 
mation, DM halos of Milky Way-type galaxies should 
contain numerous DM subhalos orbiting around the 
center, with masses between a cut-off scale 10~^^ — 
10"^ and IO^^Mq 3], where Mq denotes the solar 
mass unit. Within the resolved mass scales, numer- 
ical high resolution A^-body simulations of structure 
formation, such as the Aquarius Project [3] or the Via 
Lactea II simulation , predict the subhalos to follow 
a power-law distribution in mass, dN/dM cx M~", 
where a G [1.9; 2.0]; their density profiles resemble the 
host halo's, resulting in high central densities. Spa- 
tially, subhalos follow an "anti-biased" distribution. 



i.e., the majority orbits far away from the host's cen- 
ter. 

With lacking baryonic but dominating WIMP con- 
tent, subhalos could appear as non- variable, faint 
point-like or moderately extended high-energy (HE) 
7-ray sources without astrophysical counterpart in 
any other wavelength band. A small fraction could 
then be detectable with current high- or very-high- 
energy (VHE) 7-ray telescopes H 111 such as 
Fermi-LAT (20 MeV - 300 GeV) fiq and imaging air 
Cherenkov telescopes (lACTs; E > 100 GeV) tUj]. 

This article summarizes the main results of a study 
on the detectability of DM subhalos with Fermi-LAT , 
their properties, and the investigation of the most 
promising object from first searches for subhalo can- 
didates in the first Fermi-hAT point-source catalog 
(IFGL), IFGL J0030.7-f0724. For details and fur- 
ther explanation the reader is referred to Zechlin et 
al., 2011 [12,]. 



2. Gamma rays from DM subhalos 

The total rate of photons with energy E in the in- 
tervall [Ei; E2], originating from self-annihilating DM 
in a DM subhalo, is 



C = ^-^ IdVp (r), N, 



dx^ 
dx 



Ei/r. 



(1) 

where {<tv)cS is the thermally-averaged annihilation 
cross section times the relative velocity, the WIMP 
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mass, and dN-^/dx, x = E/m^, denotes the differen- 
tial spectrum of photons per annihilation [T3| . The 
produced photon flux is (f> = C/(At:D'^), assuming 
Ts ^ £> (see below). Consistent with numerical sim- 
ulations, the density profile of DM subhalos is taken 
to follow a Navarro- Frenk- White (NFW) profile [14], 
scaled to the characteristic inner radius rs and char- 
acteristic inner density ps- For subhalos, where (e.g., 
tidal) disturbances caused by their formation history 
are negligible, the profile parameters are related by 
the virial halo mass Afvir- The quantity Mvir is given 
by the mass inside the sphere of radius i?vir, which 
encloses a mean density of Ac times the critical den- 
sit y o f the Universe Pcrit at the considered redshift 
z [III, Afvir := 47r/3 Ac/9criti?vir- The virial over- 



density at z = is Ac ~ 100 |15|. Integrating 



the NFW profile, the subhalo mass M is given by 
Af = 47rps?'^/(c), where /(c) = In(l-hc) -c/(l-hc) and 
c denotes the concentration parameter of the subhalo. 
For non-disturbed halos, the concentration is defined 
by Cvir = Rvir/rs, depending on the subhalo mass and 
redshift. In the following, we use the low- mass ex- 
trapolation of a toy model for Cvir [13 > based on [l3] • 
Corrected by a factor depending on the galactocen- 
tric distance, effects of subhalo formation processes 
will be included. In this model, henceforth subhalo 
model (SHM), the concentration increases with de- 
creasing galactocentric distance, as indicated by nu- 
merical simulations [7] . Intrinsic to the stochastic pro- 
cess of halo formation, the concentration of individual 
subhalos scatters around the median, see [13, [Si • A 
general study of the concentration-model dependence 
and further details are given in [T2j . 

Given the relations above, Eq. [1] simplifies to 



C 



((Tw)off iV^AcPcrit Mcl 
/(Cvii 



18m2 



(2) 



Conveniently, (crw)cff is normalized to {crv)o = 3 x 
10~2^ cm'^s"^, which coincides with the correct relic 
density. A higher annihilation rate, a so-called boost 
factor ((7w)cff/(o"y)oj could in general be related to 
sub-substructures or the particular particle physics 
framework. 



3. Fermi-LfiJ sources as DM subhalos 

Based on the observational quantities flux and an- 
gular extent, candidate sources for DM subhalos will 
be selected. With Eq. [51 the effective self-annihilation 
cross section {crv)cS required for a given flux and 
intrinsic source extent 9s is determined with C — 
AirD'^cj), where 9s constrains the distance D to the 
subhalo. The intrinsic extent of a subhalo is traced 
by the characteristic profile radius Vg, since 87.5% 
of the total 7-ray luminosity is produced within rg. 



Hence, D « rs/9s, where 9s denotes the angle corre- 
sponding to . About 68% of the 7-ray luminosity is 
emitted within 6*68 — 0.46 0s- The relations below are 
given with respect to 9s, but can easily be adapted for 
9q8 , which represents the extent for comparision with 
observational data. With the WIMP model, {av)eS 
is fully determined by the subhalo mass and the ob- 
served flux and extent: 



4Ac/9c 



(3) 

For heavy WIMPs, a detection of DM subhalos is fa- 
vored in the high-energy band of Fermi-hAT (10 — 
100 GeV), due to improving sensitivity [l3| and the en- 
ergy spectrum of DM annihilation. As will be demon- 
strated later, a fiducial subhalo candidate source is 
characterized by faintness and a moderate angular 
extent. Therefore, the flux of the fiducial source 
is chosen to be at the level of detection sensitivity 
(one year), which has been accurately studied in [12]: 
(j)^'^{10 - 100 GeV) = 1.6 X lO-i^cm-^s-i for 9^'^ = 
1°, corresponding to 0|g w 0.5°; the Galactic position 
has been chosen to match IFGL J0030. 7+0724. 

For WIMPs of = 150 GcV fully annihflating 
to T+r^, the effective boost factors (crw)off /(o'w)o re- 
quired to generate the emission (j)^'^ of the fiducial 
source by DM annihilation are depicted in Fig. [Tj 
Within the halo-to-halo scatter of the concentration, 
the required boost spans about one order of magni- 
tude. Minimal boost is required for massive subhalos 
between 10^ and 10^ M©. In Fig. [TJ the necessary 
effective annihilation cross section is compared with 
current observational limits on {av)cff / {(Jv)o [Tol. [20|. 
The Figure shows that, within the scatter, the boost 
factor needed to explain 7-ray sources such as the 
fiducial with massive DM subhalos (-- 10^ - 10^ Mq) 
is consistent with observational constraints. These 
subhalos are in corresponding distances from 0.5 to 
lOkpc. Further reduction of the boost could be ac- 
complished by additional sub-substructure (factor of 
2-3, [2l|). For the considered WIMP model, this leads 
to a required boost of order unity within the scatter. 

For this WIMP model and current observational 
constraints on {av)cs, on average 0.2 subhalos be- 
tween 10^ and 10® Mq are expected for detection with 
Fermi-LAT in one year, see [l^. Given the Poisson 
distribution of the average, up to one massive subhalo 
is expected in the one-year data set. 



4. Searches for DM subhalos in the 1 FGL 

As shown above, DM subhalos would appear as 
faint, non- variable, and moderately extended sources 
without astrophysical counterparts. First candidates 
have been searched for in the 11-months point-source 
catalog of Fermi-LAT (IFGL, 100 MeV - 100 GeV, 
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Figure 1: Boost factor (o-u)cff required for the fiducial 
Fermi-LAT source to originate from a DM sublialo of 
mass M. Tiie solid green line indicates the prediction of 
the average SHM model, while its scatter is depicted by 
the dotted green lines (and the green-shaded area). 
Current observational constraints on (CTt;)off are shown by 
the light grey- and grey-shaded area. 

[23|). listing 1451 7-ray sources together with flux, 
position, significance of variability, and spectral cur- 
vature. Among the sources, 630 objects are not confi- 
dently associated with known sources at other wave- 
lengths. 

To select possible subhalo candidates, the sample 
of unassociated sources has been scanned for non- 
variable sources detected between 10 and 100 GeV. 
A high-energy detection assures subhalo candidates 
driven by heavy WIMPs and furthermore avoids con- 
fusion with high-energy pulsars [23]. The position in 
the Galaxy has been constrained to galactic latitudes 
\b\ > 20°, eliminating confusion with Galactic sources. 

Twelve unassociated sources pass these cuts. Apart 
from one, the sample consists of sources at the faint 
end of the entire IFGL sample. A further discus- 
sion of these objects is presented in [l^]. Accord- 
ing to the IFGL catalog, the counterpart search 
has been extended to a wider choice of astronomi- 
cal catalogs. Given the twelve candidates, the most 
promising has been selected, governed by lacking 
association, faintness, and spectral shape, namely 
IFGL J0030. 7-1-0724. 

Both the high-energy flux of IFGL J0030. 7+0724, 
(/)p(10 - 100 GcV) = (1.5 ± 0.7) X 10-i"cm-2s-i, 
as well as the spectral index, F = 1.7 ± 0.4, are 
well compatible with a self-annihilating DM origin. 
Note that the source has only been significantly de- 
tected between 10 and 100 GeV. Within the posi- 
tional uncertainty of the 7-ray signal (at 95% confi- 
dence), IFGL J0030.7-f0724 has no convincing muhi- 
wavelength counterpart. Furthermore, X-ray observa- 
tions with ROSAT [IJ] constrain the energy-fiux of 
a possible X-ray counterpart to < 10^^^ crgcm^^ s^^ 



4.1 . Fermi-LAT data 

For further investigation of possible counterparts, 
temporal variability, and the angular extent of 
IFGL J0030.7+0724, updated results based on the 
24-months public archival data set of Fermi-LAT 
between 10 and 100 GeV have been derived. The 
point-source analysis has been performed with the 
Fermi Science Tools v9rl8p6 [26j . using recom- 
mended options and the instrument-response func- 
tions P6_V3.DIFFUSE Details are given in [li|. 

Within a radius of 0.5° around the nominal source 
position, six photons (one class 3 event, five class 
4 events) between 10 and 100 GeV have been de- 
tected after 24 months. With respect to the signal, 
the influence of the Galactic foreground as well as 
the extragalactic background at the source position 
is negligible. The integrated flux between 10 and 
100 GeV, reconstructed with the point-source analy- 
sis, is (/)p(10-100GeV) = (0.9±0.4) x 10-^° cm-^ g-i. 
The significance of the signal in this energy-bin is 
about 6.6 Gaussian standard deviations. 

The variability of IFGL J0030. 7-1-0724 has been 
tested for, analyzing the temporal photon distribu- 
tion for compatibility with a constant flux with an un- 
binned Kolmogorov-Smirnov test. The test confirms 
the null-hypothesis of a steady flux with a probabil- 
ity of about 0.5. Given the positional distribution of 
the signal photons, the (intrinsic) spatial extent has 
been tested with a likelihood-ratio test, comparing the 
null-hypothesis of a point-source with the assumption 
of the intensity distributed following the line-of-sight 
integral of the squared NFW proflle. The data favors 
a moderate extent 9s = 0.14^q;^2 deg, which is, how- 
ever, not signiflcantly different from a point-source hy- 
pothesis. The upper limit is 6*5 < 0.72 deg at 95% 
confidence level. 

4.2. Sw/ff-XRT data 

New observations of the field with the X-ray tele- 
scopefXRT, 0.2-10 keV) onboard the Swift sateL 
lite [2g] with a total effective exposure of 10.1 ks 
led to the discovery of seven new X-ray sources. 
The (unabsorbed) energy fluxes of the discovered 
sources range from ~ 2 x 10~^^ ergcm"^ s~^ to 2 x 
10-13 erg cm-^s-i between 0.2 and 2keV. See for 
further details. 



5. Discussion 

5.1. IFGL J0030.7+0724 as AGN 

Within the updated positional information on 
IFGL J0030. 7-1-0724 (see Sect. IH]), the radio 
source NVSS J003119-H072456 (/i.4GHz = (11-6 ± 
0.6) mjy) provides a viable counterpart. The radio 
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source positionally coincides with the newly discov- 
ered hard X-ray source SWIFT J003119.8-h072454 
(r = 1.6 ± 0.3), energy flux ~ 2 x lO'^^ gj.gcm-2 s'^ 
between 0.2 and 2 keV, and the optical counterpart 
SDSS J003119.71-h072453.5 (m(r) = 17.4'"). Nor- 
malized to the radio flux, the average spectral energy 
distribution of a high-energy-peaked blazar (HBL, see 
[295) fits the data, assuming standard temporal vari- 
ability of the source. Reinforcing, the HBL scenario 
is favored by the spectral indices of the X- and 7-ray 
signal. 

The optical counterpart appears point-like with 
m{R) = 18.6™ (the SDSS source coincides with USNO 
0974-0005617). Under the assumption of an elliptical 
host galaxy with Mn = —23.1™, typical for blazars 
[30| . the corresponding distance modulus implies a 
redshift of z = 0.39 ± 0.03 (K-correction). 

5.2. 1FGL J0030.7+0724 as DM subhalo 

Given the absence of clear indication for variability, 
a DM subhalo origin of IFGL J0030. 7+0724 remains 
plausible. The flux and the upper limits on the an- 
gular extent provided in Sect. 14.11 are consistent with 
a subhalo of mass between 10^ and 10^ Mq. Assuming 
a mass of 10^ M0, the distance would be 2.4j:J;?kpc 
given the scatter of the concentration model, see Sects. 
H and H In case of a WIMP of 150 GeV annihilating 
to r+T^, the required minimal boost is 3 for a high- 
concentrated subhalo with a corresponding distance of 
1.7 kpc, and 13 for an average-concentrated subhalo 
with a corresponding distance of 2.4 kpc. An even 
lower boost may be required, if the subhalo contains 
sub-substructures or exhibits a cuspier profile. 



6. Summary and Outlook 

Wc have presented a state-of-the-art investigation 
on the detectability of DM subhalos with the Fermi- 
LAT. Concluding, one massive DM subhalo will prob- 
ably appear in the first catalog. First searches have 
indeed revealed twelve possible candidates, where the 
most promising, IFGL J0030. 7-1-0724, has been stud- 
ied in detail. A detection of temporal variability and 
improved astrometry would clarify the physical origin 
of this object. 

The accompanying paper Zechlin et al., 2011 JJiJ 
will contain additional details. The recent release of 
the second Fermi-LAT catalog, 2FGL, allows a deeper 
search for subhalo candidates, which will be presented 
in a future work. Note that related work is presented 
in [3I] and [ssj . 
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